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1. Introduction: Reasons for physics beyond the Standard Model 

Although the Standard Model (SM) of particle physics provides an excellent description of 
electroweak and strong interactions, there are many reasons that we expect to observe new forces 
giving rise to new particles at larger masses than the known fermions or bosons. One oft noted 
source of this belief is the observation of dark matter in the cosmos as evidenced by galactic angular 
velocity distributions [1], gravitational lensing [2], and galactic collisions [3]. The existence of dark 
energy, believed to cause the accelerating expansion of the Universe, is another source of mystery 
[4]. The fine tuning of quantum corrections needed to keep, for example, the Higgs boson mass at 
the electroweak scale rather than near the Planck scale is another reason habitually mentioned for 
new physics (NP) and is usually called "the hierarchy problem" [5]. 

It is interesting to note that the above cited reasons are all tied in one way or another to 
gravity. Dai^k matter may or may not have purely gravitational interactions, dark energy may be 
explained by a cosmological constant or at least be a purely general relativistic phenomena, and the 
Planck scale is defined by gravity; other scales may exist at much lower energies, so the quantum 
corrections could be much smaller. There are, however, many observations that are not explained 
by the SM, and have nothing to do with gravity, as far as we know. Consider the size of the quark 
mixing matrix (CKM) elements [6] and also the neutrino mixing matrix (PMNS) elements [7]. 
These are shown pictorially in Fig. 1 . We do not understand the relative sizes of these values or nor 
the relationship between quarks and neutrinos. 
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Figure 1: (left) Sizes of the the CKM matrix elements for quark mixing, and (right) the PMNS matrix 
elements for neutrino mixing. The area of the squares represents the square of the matrix elements. 



We also do not understand the masses of the fundamental matter constituents, the quarks and 
leptons. Not only are they not predicted, but also the relationships among them are not understood. 
These masses, shown in Fig. 2, span 12 orders of magnitude [7]. There may be a connections 
between the mass values and the values of the mixing matrix elements, but thus far no connection 
besides simple numerology exists. 

What we are seeking is a new theoretical explanation of the above mentioned facts. Of course, 
any new model must explain all the data, so that any one measurement could confound a model. 
It is not a good plan, however, to try and find only one discrepancy; experiment must determine a 
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Figure 2: Lepton and quark masses. 



consistent pattern of deviations to restrict possible theoretical explanations, and to be sure the new 
phenomena is real. 

2. Use of flavour physics as a new physics discovery tool 

While measurements of CKM parameters and masses are important, as they are fundamental 
constants of nature, the main purpose of flavour physics is to find and/or define the properties of 
physics beyond the SM. Flavour physics probes large mass scales via virtual quantum loops. Ex- 
amples of the importance of such loops are changes in the W boson mass {M\^) from the existence 
of the t quark of mass rut, dMw oc rnj, and changes due to the existence of the Higgs boson of mass 
Mh, dMw °^ ln(M//). 

Strong constraints on NP are provided by individual processes. Each process provides a dif- 
ferent constraint. Consider for example the inclusive decay b — sj. The SM diagrams are shown 
in Fig. 3(a). The SM calculation considers either a B or meson and then sums over all decays 
where a hard photon emerges. The Feynman diagrams corresponding to a NP process with a virtual 
charged Higgs boson are shown in Fig. 3(b). 




Figure 3: (a) SM diagrams for the quark level process b ^ sy and (b) NP diagrams mediated by a charged 
Higgs boson. 
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My average of measured branching fractions, including a recent result at this conference from 
BaBar [8], is 

^(Z7 ^ 57) = (3.37 ±0.23) X 10"^ (2.1) 

The SM prediction is (3.15 ±0.23) x 10""^ [9, 10]; thus the ratio of the measured to theoretical 
prediction is 1.07±0.10, which serves to severely limit many, if not most, NP models. An exam- 
ple of such a constraint for the two-Higgs-double model (2HDM), with the ratio of the vacuum 
expectation values of the two Higgs doublets having a value of tan j3 = 2, is given in Fig. 4. The 
intersection of the lower edge of the NP prediction with the upper line of the experiment gives a 
limit at the approximately two standard deviation (a) level of m//+ > 385 GeV [1 1]. 
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Figure 4: .^{b sy) as a function of the charged Higgs boson mass. There are three sets of curves including 
the experimental average shown in dotted (red), the SM prediction in dashed (blue) and the 2HDM model in 
solid (black). The middle curve of each set gives the central value and higher and lower ones the ±1(7 bands 
[9]. 

It is of course true that experiments measure the sum of SM and NP, so to set limits or see 
signals from NP we have to be confident in understanding the SM contribution. This was clear in 
the above example of — )• sy. In other processes SM calculations are not possible, or are not un- 
ambiguous. An often used paradigm is that "tree-level" diagrams are dominated by SM processes, 
but loop level diagrams can be strongly influenced by NP. Examples of tree diagrams are given in 
Fig. 5(a), while both SM and loop diagrams that could have new as let undiscovered particles in the 
loop are shown in Fig. 5(b); here both diagrams would add in amplitude and thus interfere allowing 
for either increases or decrease in rates and asymmetries. 




a) SM tree diagrams b) SM loop diagram c) NP loop diagram 



Figure 5: (a) Examples of SM tree-level diagrams, (b) Example of a SM loop diagram for mixing 
interfering that can interfere with (c) a NP diagram. 
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NP can also be limited via a general study of loop processes, again assuming that it doesn't 
contribute to tree-level processes. One can write a general local Lagrangian which includes an 
infinite tower of operators with dimension d > A, constructed in terms of SM fields, suppressed by 
inverse powers of an effective scale A > M-w as 

(d) 

^eff = ifsM+I^^i'\ (2.2) 

where the ^'s represent the SM fields [12]. Lower bounds on A are shown in Table 1, taking the c/ 
coefficients equal to one for different operators. 



Operator 


A lower bound (TeV) 


Observable 




9.8 X 10^ 


Kl — Ks mass difference, AniK 


\m\{sLtdLY] 


1.6 X 10"^ 


CP violation in Kl decay, Ek 


R&[{sRdL){sLdR)\ 


1.8 X lO'^ 


Kl — Ks mass difference, AniK 


lm[{sRdL){sLdR)] 


3.2 X 10^ 


CP violation in Kl decay, Ek 


MibttdLY] 


9.8 X 102 


Bl — Bs mass difference, Amg 


ImlibLrdtf] 


1.6 X 10^ 


CP violation in — ;> J/^Ks decay 


^e\{bRdL){-SLdR)] 


1.8 X 10^ 


Bl — Bs mass difference, Amg 


\m[{bRdL){sLdR)] 


3.2 X 10^ 


CP violation in B" — ^ J/^Ks decay 



Table 1: Lower limits on the size of representative dimension six loop diagram operators (A), taking the 
coefficients c, equal to one. Adopted from Ref. [12]. 

These limits are impressive ranging from approximately 100 to 10,000 TeV, well beyond the 
direct collider limit search possibilities. Of course it is possible that the c, are smaller than one, 
but they have to be in range from 10^^ to 10^^ to bring these limits to the TeV scale. There are 
ways, however, that NP can be consistent with these limits. One way is for the new particles in 
the loop diagrams to be degenerate in mass and have opposite phase, so they cancel in the loops. 
Another is for the quark mixing angles to be the same in the NP sector as in the SM. This defacto 
prescription, called "Minimum Flavor Violation," is well described by Buras [13]. Of course very 
heavy new particles are also a possibility. Overall, these measurements put very severe constraints 
on NP models. 

3. IMixing and CP violation 

Let us consider in more detail the mixing diagrams which allow transformation of neutral 
mesons into their anti-particles. Fig. 6 shows the second order weak interaction diagrams for D*', 
S° and S? mesons. The amplitude strength depends on the values of the CKM elements at the four 
vertices and approximately as the magnitude of the mass-squared of the virtual quarks permitted. 
B" and B° mixing transitions are far larger than that of the D", because the t quark is allowed in 
the B cases, but for the D° the heaviest allowed exchange quark is the b. Hence D" mixing is quite 
small, but finite [28], while the B*' is quite significant and the B'' huge, the difference in the two B 
cases due entirely to relative sizes of the CKM elements. 

The Schrodinger equation is used to describe time dependent transformations between meson 
and anti-meson states for diagrams similar to the one in Fig. 5(b) in terms of a mass matrix M and 
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a decay matrix T: 




Mn-/rn/2Mi2- 

M*2-/r*2/2 M22- 




(3.1) 



Diagonalizing we find \Ml) = p \B^) + q \B^), \Mh) = p — q where p and ^7 are complex 
numbers that are equal to 1 / \/2 in the absence of CP violation. Also, assuming CPT invariance, the 
mass m{Bs) = Mn = M22 = {Mh +Mi) /2, and the mass difference is defined as AM = Mh—Ml, 
while for the decay widths, T = 1 /t(B,) = [Th + Tl) /2 = Tn = r22, and AF = Tl - T//. For 
convenience the variable y = Ar/2r is defined. 

Time dependent CP violation depends on the process in many ways. If the final state / is a CP 
eigenstate then f = f and the CP violating asymmetry is defined as [16] 

r(M^/)-r(M^/) 

alt) = — 7= 7 ; . (3.2) 

r(M^/)+r(M^/) 

The amplitudes are defined as Ay = A(M—)'/),Aj=A(M —)•/). The magnitude of the CP violating 
effect for each channel depends on the CKM elements present in the decay, and is given by 

A, = ^^. (3.3) 

qAf 

The decay rates, when = 1, are given by 

r(M^/) = Nf\Af\^e^^' ^cosh^-ReA/sinh^+ImA/sin(AM?; 

cosh — -ReA/sinh—- ImA/ sin (AMOj . (3.4) 

Thus the CP violating asymmetry is given by 

ImAfsin(AMO 

ait) = tt/ ^ —KrT ■ (3-5) 

cosh^-ReA/sinh^ 
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3.1 CP violation in decays 

Let us consider the decays — )• J/y^^ where (j) — )• K^K and J/xf/^^^ - The tree-level 
decay diagram is shown in Fig. 7. Here we use the interference between the direct decay and the 
one that proceeds via mixing: B° — — ^ J/y^^ for example. The CP violating asymmetry a{t) 




Figure 7: Tree level diagram for B^. -^• J/\jf7t^7t or K^K decays. 



can be expressed in terms of the phase angle (^s, and the CP parity of the final state, T]cp- In the SM 
Im(A/) = -77cpsin(^>, = Tjcparg (^^^) = -0.04 rad for CP odd states [17]. 



V,1,V,i,Vc,V*J 

Based on a conjecture by Stone and Zhang [18], the LHCb collaboration found the decay 
mode B^. — > 7/va/o(980) where /o(980) K^K^, close to the predicted level [19], and was soon 
confirmed by others [20]. As the final state consists of a scalar /o(980) and a vector J/y, the 
relative angular momentum between the /o(980) and the is a P-wave and the final state is pure 
CP odd. LHCb then examined the entire tt+tt^ mass spectrum in the 5^ — J/\j/K^7i^ channel, 
shown in Fig. 8. The region between the arrows has a large signal to background ratio. In this 
region a full Dalitiz like analysis was performed that identified the individual components, and 
more importantly showed that the final state is > 97.7% CP-odd [21]. Measurement of CP violation 



based on Eq. 3.5 gives a value 0, = —0.019 



+0.173+0.004 
-0.174-0 



oo5 rad [22]. 




500 



1000 1500 2000 

m(jt+Jt") (MeV) 



Figure 8: The n^n invariant mass in — > J/xittz^tz decays. The dashed (red) curve shows the back- 
ground, and the arrows indicated the region chosen for CP studies. 
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The J/y^ state being composed of two vector particles is not a CP eigenstate. Such final 
states still can be used but an angular analysis is necessary to separate the CP-even and CP-odd 
components [23], which requires fitting time dependent angular distributions. The three angles are 
shown in Fig. 9. 




Figure 9: Pictoral description of the decay angles. On the left 6 and defined in the J/yf rest frame and on 
the right \j/ defined in the ([) rest frame. (From T. Kuhr [24].). 



In addition to the three P-wave amplitudes in the decay width, an additional amplitude due 
to possible K^K S-wave also must be considered [18]. The complexity of the analysis scheme 
is shown by writing the distribution of the signal decay time and angles, described by a sum of 
ten terms, corresponding to the four polarization amplitudes and their interference terms. Each of 
these is the product of a time-dependent function and an angular function [23] 

oc y ht(t)fk(e,w,d)). (3.6) 

dtded\i/d<p 

The time-dependent functions hj^{t) can be written as 

hkit) = Nte^^"' [c<:Cos(Ani.v?) + Ji.sin(Anii?) +a^cosh (jAr^f) +Zji.sinh (iAFv?)] . (3.7) 

where Am., is the B° oscillation frequency. The coefficients A'jt and ak,---,dk can be expressed 
in terms of 0v and four complex transversity amplitudes Aj at ? = 0. The label / takes the values 
{_L, 11,0} for the three P-wave amplitudes and S for the S-wave amplitude. For a particle produced 
in a flavour eigenstate the coefficients in Eq. 3.7 and the angular functions fi^{6,\l/,<p) are then, 
given in Table 3.1 [26], where 5o = is chosen arbitrarily 

The differential decay rates for a meson produced at time ? = are obtained by changing the 
signs of ips, ^_l(0) and As (0), or, equivalently, the signs of and dt-^ Note, that for the J /'^ifK^K^ 
final state, only line #7 in Table 3.1 contributes. 

Results of the analyses are presented in the AF^ — <ps plane, since earlier low statistics data had 
significant coiTclations. The current results are shown in Fig. 10 for several experiments. Th LHCb 

'The decay width is invariant under the transformation (0j, AF^, 5|| , 5x, 5s) t-^ {n — (^j, — AFj, — 5|| , ;r — 8^,— 5s) 
which gives rise to a two-fold ambiguity in the results. This ambiguity has been removed by analyzing the interference 
between the S and P waves [27]. 
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Table 2: The components of the decay width contributing to Eqs. 3.6 and 3.7, where 5,^ = 5, — 5j, e.g. 

5x0 5x - 5o. 



result, 0s = 0.001 ±0.101 ±0.027 rad is the most accurate and dominates the average. Combining 
with the J/\ifK+n^ result LHCb finds [25] 



0^ = -0.002 ± 0.083 ± 0.027 rad, 
which is consistent with the SM prediction. 

LHCb 1 .0 fb"' + CDF 9.6 fb + D0 8 fb"' + ATLAS 4.9 fb"' 
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Figure 10: Measurements of AT versus 0^ as compiled by HFAG [28] 



I choose to determine the lifetime and AT, based on measurements of only fully recon- 
structed decays where all the final state particles are seen. Fig. 1 1 shows the data from four separate 
determinations using //v^0, two using //i/a/o(980), one using D^n^ at fully mixed CP state, and 
one using K^K .The. last mode is a CP-even eigenstate where the effects of CP violation are taken 
to be negligible. The results to an overall fit are shown in the black oval covering 39% confidence 
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level (CL). They are 

= 1.509 ± 0.010 ps, 

0.092 ± 0.01 1 ps"^ 
AT, 

These result agree very well with the theoretical predictions. 



0.07 ±0.01 



(3.9) 




1.3 1.4 1.5 1.6 1.7 

l/Fs [ps] 



Figure 11: Measurements of ATi versus the effective hfetime in each mode labeled as 1 /F^, and the average 
as shown in the black ellipse as determined by A. Phan [29]. The ovals show 39% CL contours, while the 
bands are at 68% CL. 



There is another way CP violation can show itself in decays. Consider the asymmetry 

defined as _ 

Y{M^f)-T{M^f) 



ait) 



T[M^f)+T{M^f) 



(3.10) 



The difference here with Eq. 3.2 is that the final state / is flavour specific, / 7^ /, and at zero 
decay time both T [M —;•/)= and T {M — )• /) =0. For example, in the case of decays the 
asymmetry between the decay rates forS^ -^D^jx^v andB° — )• Dj }X^v can be measured. Another 
method is to use events where both Z^-flavoured hadrons decay semileptonically, i.e. toXpLV. When 
one b mixes and the other decays these events produce like-sign muons, and the difference between 
jU+/i+ and iti^iU^ events can be examined. 

In the SM the asymmetry is related to decay width as fly/ = (Ar/AM)tan^i2, where tan ^12 = 
— arg(— ri2/Mi2). The asymmetries are expected to be very small in the SM, -4.1 x 10""^, and 
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1.9 X 10^^ for BP (a'^i) and («J;)> respectively [30]. The DO collaboration, however, reported 
an anomalously lai^ge value for the dimuon asymmetry of A^^ = (—0.787 it 0.172 it 0.093)% [31]. 
As they are summing over and decays, they graph their results as the purple band shown in 
Fig. 12. As the width of the band corresponds to itla uncertainty, their result is 3.9a from the SM 
value. If confirmed, this would be a clear indication of NP. 




-0.04 -0.02 0.02 



Figure 12: Measurements of the semileptonic asymmetry using dimuon events from DO. The purple band 
(A'^i) is their base result, while the bands have impact parameter (IP) cuts applied that change the relative 
andfiO yields [32]. 



Subsequently DO measured the individual semleptonic asymmetries a^i [33] and a^^ [34] using 
semileptonic B decays, where a single muon is detected in conjunction with a charmed meson. 
The invariant K^K^n"^ mass distribution, where the K^K is required to be consistent with the (j) 
meson mass, for putative B^ -^Dfpi^X decays is shown in Fig. 13, where theX indicates a missing 
neutrino plus possibly additional mesons. The events in the Ds signal peak come predominantly 
from B° decays. The measurement results from counting the number of Dfn versus events 
and correcting for any residual detector asymmetries. From these data DO finds a^j = (— 1.08 ± 
0.72 ±0.17)%, and from studying D^/i^x decays, with D± K^n+n-, af, = (0.68 ± 0.45 ± 
0.14)%. 

Combining their di-muon results, including some separation into a^f and a^^ values based on 
muon impact parameter, with the semileptonic measurements, Bertram compares the DO measure- 
ments with the SM in Fig. 14 [35]. 

There is a differing world view on these asymmetries. The B -factories have produced mea- 
surements of a'lj that average to 0.0002 it 0.0031 [28]. LHCb recently measured the asymmetry in 
semileptonic B° decays, in a similar manner as DO, using a 1.0 fb"' data sample. The K+K n"^ 
mass spectra is shown in Fig. 15 for candidate decays with a muon of opposite charge to the can- 
didate Ds meson. LHCb periodically reverses the polarity of their spectrometer magnet, as did DO. 
This data is for one configuration only (magnetic field down). The excellent mass resolution al- 
lows the D+ and D+ signals to be completely separated. LHCb has different detection efficiencies. 
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Figure 13: The mass spectrum for events where the tracks form a common vertex detached from the 
primary. The peaks are at the and D+ masses. (Note the zero suppressed vertical scale.) 




Figure 14: Combination DO semileptonic decay asymmetry measurements with measurements using di- 
muons separated by into two impact parameter categories by a cut of 120 /zm. The error bands represent Id 
uncertainties on each individual measurement. The ellipses represent the 1, 2, 3, and 4c7 two-dimensional 
C.L. regions, respectively. The SM point is shown with a black dot. 



in principle, for ;r+ versus and pL^ versus jx . In each magnet setting the tracking efficiency 
differences mostly cancel considering the binary pairs n^pL^ versus jx^ . Most of the remain- 
ing differences are cancelled by averaging the two opposite magnet poIarities.The different muon 
triggering efficiencies are evaluated using a sample of about one million jU+/i^ found in 

events that were triggered independently of the J /^. LHCb finds = (—0.24 ±0.54 it 0.33)% 
[36]. Measurements from the B -factories and LHCb are plotted in Fig. 16. The data are completely 
consistent with the SM, and differ somewhat from the DO results. 
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Figure 15: Invariant mass distributions for (a) K^K 7r+ candidates and (b) K^K^ candidates for magnet 
down data, requiring that m{K^K^) be within ±20 MeV of the meson mass. The vertical dashed lines 
indicate the signal region. 
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Figure 16: The LHCb measurement of a'^, and the B-factory measurement of a^^. The band edges are the 
central values ±1(7 uncertainties. 
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3.2 Mixing and CP violation in Charm 

LHCb recently measured mixing in the — T^ system at more than 5 a significance [37], 
consistent with the average of several other less significant results [28]. Measurement of charm 
mixing at the levels seen is expected in the SM. Measurement of CP violation at the ~ 1% level, 
however, was not expected. An interesting paper discussed the possibility of seeing NP in charm 
CP violation specifically by looking at singly Cabibbo suppressed decays [38]. 

Using the same definition for asymmetry as in Eq. 3.2, the Cabibbo suppressed final states 
D" — K^K and — > tt+tt^ have been investigated. A flavor tag is provided by using D" that 
result from D** — )• tt^D" decays. In order to cancel the effects of D** production asymmetries and 
71+ versus detection asymmetries, the difference between time-integrated asymmetries AAcp = 
acp{K^ K^) — acp{'Jl^ n^) is determined. LHCb first published a 3.5 standard deviation result [39]. 
When averaging data from CDF [40] and Belle [41], 1 find AAcp = (-0.74 ±0.15)%. However, 
none of the experiments by itself has a result of more than 3.5a significance, and more data are 
needed to see if the result is real. Meanwhile the theory community continues to argue if this 
asymmetry difference can be explained by the SM [42] . 1 choose to treat this as a limit on NP of 
1%> -AAcp > 0%. 

4. Exclusive rare decays 

We have seen that the inclusive "rare" decay b ^ sy has a relatively accurate SM prediction 
(see Sec. 2). Other well determined predictions exist for exclusive rare decays. Since these de- 
cays have suppressed branching fractions the possibility that NP can make large contributions is 
enhanced. 

4.1 

Heavy suppressed in the SM partially due to helicity conservation, the decay — jLt+/i^ 
is predicted by Buras [43] to have a branching fraction, after a correction due to the relatively 
large size of AF^, of (3.5 ±0.2) x 10^^ [44]. The SM branching fraction for the corresponding 
b" decay is even smaller due to the smallness of the CKM element \Vtd\ compared with \ Vts\- The 
SM diagrams are shown in Fig. 17, along with possible interfering diagrams from NP processes. 
We note that many NP models are possible not just the one example of super-symmetry that is 
shown. While upper limits on this process have been established by several experiments, as shown 
in Fig. 18, a new LHCb result sees evidence of a signal [45].^ 

LHCb has unique capabilities at the LHC to study purely hadronic decays of Z^-flavoured 
hadrons. In the case of this analysis such decays are used to study selection criteria, backgrounds 
and to measure the transverse momentum, pT, dependence of the ratio of to B° production, 
called fslfd- Multivariate discriminants are trained using several variables, including track impact 
parameters, B lifetime, B pr, B isolation, muon isolation, minimum impact parameter of muons, 
etc... The signal properties are established using fully reconstructed B decays into two oppositely 
charged hadrons. The signals for four different final states are shown in Fig. 19. 

^Most LHCb results reviewed here use a 1.0 fb^' data sample. This result uses 2.1 fb^' . 
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Figure 18: Upper limits for — > from different experiments, as of July 4, 2012. The "LHC combi- 

nation," includes individual ATLAS, CMS and LHCb contributions shown in the figure. 



The analysis is done applying two different multivariate discriminates in sequence. The last 
one is a Boosted Decision Tree (BDT), which is constructed to have the signal probability flat in 
BDT, and the background peaked around zero. The data show a 3.5 standard deviation significance 
signal from the final fit. The projection for the final boosted decision tree variable, BDT>0.7 is 
shown in Fig. 20. 

To find a branching fraction one needs to determine not only the number of signal events and 
their average detection efficiency, but also the total number of events. This number is inferred by 
measuring the number of events in the similar decay channels B — )• J/\j/K^ and — )• K^Ti^, us- 
ing their known branching fractions, and knowledge of the production ratio for B^JB^ decay, fs/ fd 
(assuming that B and B° production are equal.) LHCb determines this ratio from two sources. 
The most precise is semileptonic b decays [46], and the second is specific hadronic channels us- 
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[MeV/c2] 

Figure 20: Invariant fi^l-i^ mass distribution of the selected candidates (black dots) with BDT > 0.7. The 
result of the fit is overlaid (blue solid line) and the different components detailed: (red long dashed), B 
(green medium dashed), h^h'^ (pink dotted), if — > n^pL^v (black short dashed) and B — > nn^jJ.^ 

(light blue dot dashed), and the combinatorial background (blue medium dashed). 



ing a theoretical model. The average of the two methods gives fs/fd = 0.256 it 0.020 [47]? The 
hadronic channels also allows an exploration of the pseudo-rapidity, rj, and pj dependence. It 
is found to be independent of rj, but in px the B^/b" yield decreases as shown in Fig. 21. For 
LHCb the px distribution of the signal is quite similar so little systematic uncertainty is introduced. 
Other LHC experiments have in the past used the LHCb value for fs/fd', they will have to take into 
account this pj dependency in the future. 

The most important implication of this measurement is that it is consistent with the SM predic- 
tion. It rules out many SUSY models especially those with large values of the parameter tan j3 [48]. 
Many NP models are eliminated or severally constrained when looked combined with other mea- 

The hadronic channels used are B D+K-, B D+K- and B, D+K'. 
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Figure 21: The variation of fs/ fd as a function of pf measured by LHCb using fully reconstructed hadronic 
decays. 



surements. For example, the concurrently available upper limit on — )• pL^pL^) < 0.94 x 10^^ 
at 95% CL is plotted versus ^{B^^ — )• }X^}X^) for several models including the SM in Fig. 22(a). 
Here the space taken by any particular model is caused by a "reasonable" variation of its param- 
eters. In Fig. 22(b) the models are plotted in versus ^(B^ — )• plane. Many models ai^e 
gone or mostly eliminated. SM4, a fourth-generation model, is eliminated by the Higgs discovery. 
Here the 4th generation quarks would cause the Higgs production cross-section to be nine times 
larger and suppress the decay into yy [49]. Other models with multiple Higgs doublets are allowed 
for heavy fourth-generation quarks with masses larger than 400 GeV [50]. 




10 20 30 40 50 -1.0 -0.5 0.0 0.5 1.0 



Figure 22: (a) Correlation between .^(B^ — > pi^li^) and ."^(if — > fi^li^) in minimum flavor violation 
(MFV) model, the four-generation SM4 model, and four SUSY models. Correlation between ^{B^ — > 
/i+jU^) and the mixing-induced CP asymmetry (j)^ in SM4, the two-Higgs doublet model with flavour blind 
phases, and three SUSY models. In both cases, the SM point is marked by a star. The gray area is ruled out 
experimentally. In both (a) and (b) the blue dashed region outlines the allowed region within one standard 
deviation of their measured values. Adapted from D. Straub [51]. 



4.2 B K*i+e- 

The process B — )• K*£'^£^ is well worth studying as it proceeds via two SM model loop di- 
agrams, thus allowing for substantial NP contributions, and it is well predicted in the SM [52]. 
The diagram for the specific decay channel B ^ K pi^pi^ is shown in Fig. 23. Although this 



17 



New physics from flavour 



Sheldon Stone 



is an exclusive hadronic decay as contrasted with corresponding inclusive decay b — )• si'^l^ , there 
is sufficient precision in the theoretical predictions to make this an exceeding interestly mode to 
study [52], especially useful as the exclusive decay is much more experimentally accessible. In- 
deed, there are many suggestions for variables to use to uncover NP in this decay. Here 1 will use 
just the four- momentum transfer between the B and the K* , and the forward-backward asymme- 
try, Afb, defined as the angle of the /i+ with respect to the B direction in the di-muon rest frame. 
An essentially identical diagram can be drawn for the KpL^ jjL^ , but the spin-0 K provides a less 
interesting final state in terms of decay amplitudes. 




Figure 23: Decay diagram for B K* jX^jX . The isospin related diagram for B decays is almost 
identical, as is the one for e+e^ in the final state. 

The decay rates as a function of from four experiments BaBar [53], Belle [54], CDF [55]. 
and LHCb [56] are shown in Fig. 24 and compared with a theoretical prediction [57]. All four 
experiments agree with the prediction; LHCb has the most precise measurements. 
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Figure 24: Measured decay rates for 
Bobethef aZ. [57]. 



15 20 

[GeV^] 

p {b^ K*^IJ.^H^^ compared to a theoretical calculation from 



ApB is more sensitive to NP. The LHCb results and an average of the data from all four experi- 
ments shown in Fig. 25, however, are quite consistent with the SM prediction. LHCb also measure 
the point at which ApB crosses zero as (4.9^} 3) GeV^, where the SM predictions range from 4.0 to 
4.3GeV2 [58]. 
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Figure 25: Data and theory comparison of the forward-backward asymmetry, Afb for B —y K H^ji^ 
decay. The data are from LHCb, and the theory curves are from Bobeth et al. [57]. The SM prediction is 
shown shaded with the band covering ±1 standard deviation uncertainty estimates. The non-SM curves are 
labeled with C,, the meaning of which will be discussed in the next section. Also shown is the data averaged 
over all four experiments, where the horizontal axis is shifted somewhat with respect to the LHCb data. 



4.3 Other rare b s processes 

Other processes dominated by the b ^ s transition can be combined with the previously dis- 
cussed ones in a common analysis framework. The formalism developed starts with writing the 
effective Hamiltonian as [59] 




(4.1) 



where the C, and O, represent SM coefficients and operators, and the C- and 0\ represent NP 
coefficients and operators. The operators most sensitive to NP effects are 



O7 = '^{sa^,PRb)F^\ O, = ^-!^(s-a,yT"PRb)G^'^'', 

e 

O9 = {sy^PLb) {If I) , Oio = {sy^PLb) {Ifysf) , 

Os = m,, {sPRb) {U) , Op = mt, {sPpb) (£75^) , (4-2) 

where mfc denotes the running ^ quark mass in the IVIS scheme, Pi = {l — y^) /2, andP^? = (1 + 75)/2. 
The NP operators O' are found by switching Pj^ for Pi and vice- versa [60]. 

There are other rare b ^ s processes that contribute to the NP search. One important one is the 
exclusive decay B° — K*^y, K*^ — > Ks'nf . Time dependent CP violation for this process is given 
by 

r {f{t) ^ K*^y\ - r {B^{t) ^ K*^y) 

— ^— 1- = 5/f.ysin {^mdt) - C/f.^cos (Amj?) . (4.3) 

T[W{t)^K*^y] +Y{B^{t)^K*^y) 
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Including generic NP operators, the expected asymmetry is parameterized as 

Sk'j = , ^ , A m (e-^'^CiC'j) . (4.4) 

The SM prediction for Sk*j, where Cy = 0, is (—0.023 it 0.016), while the measurement from 
BaBar and Belle [61] is —0.16 ±0.22. Cleai^ly reducing the error would be useful, but the mea- 
surement, as we will see, is already useful. 

Another useful process is the inclusive decay b — )• sl'^t^ . Here two four-momentum transfer, 
q^, regions are defined, one at low values 1 < < 6 GeV^ which again probes C7 and the other 
with > 14.4 GeV^ which examines both C9 and Cjq. This is similar to the coefficients probed in 
the exclusive — )• K*'^l^l^ channel. 

Putting all of the channels together Altmannshofer and Straub provide constraints on the real 
and imaginary parts of the NP parameters shown in Fig. 26. The SM points are at (0,0), consistent 
with the analysis within 2a in all cases. 




RetC,) Re(C,) Re(Cjo) 

Figure 26: Individual two standard deviation constraints on the primed Wilson coefficients as well as com- 
bined one and two standard deviation constraints shown in red (dark and light) from B (brown), 
^[B X,Y) (yellow), Acp{b -J> sj) (orange), B -> /:*7(puiple), B K*n+n- (green), B K^+^- (blue) 
and Bs M^M" (gray) from [60] . 



5. Tree level decays with "abnormal" behavior 
5.1 B T V 

Here we depart from the paradigm of considering NP via loop diagrams and look instead at a 
the suppressed diagram shown in Fig. 27. The branching fraction in the SM is given by 

^{B- ^ T-v) = (1 _ ^plK^W (5.1) 

where Vub is a small CKM element and /g is the B meson decay constant. Helicity suppression, 
caused by having the spin-0 B decay into a left-handed lepton and a right-handed anti-lepton 
whose helicities align for massless leptons and thus tend to form a spin-1 system, is partially 
broken due to the relative large value of mt with respect to mg- . New physics could enter here 
by having another charged boson, e.g. a virtual H also participate in the decay instead of the 
[69]. 
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Figure 27: Feynman diagram for ^(B — > t v). 

The measurements of the branching fraction have been made by BaBar and Belle. Both col- 
laborations have separate determinations where they either fully reconstruct a hadronic decay or a 
semi-leptonic decay tag, and look for evidence for a decay on the other side requiring little extra 
energy in the event after accounting for the tag and the decay products. The experimental mea- 
surements have been compiled by Rosner and Stone [64] and listed in Table 5.1. The Belle average 



Experiment 


Tag 


^ (units of 10"^) 


Belle [65] 


Hadronic 


1 70+0.56+0.46 
^ - '^-0.49-0.51 


Belle [66] 


Semileptonic 


1 +0.38+0.29 
^■■'^-0.37-0.31 


Belle 


Our average 


1.62 ±0.40 


BaBar [67] 


Hadronic 


1.80l|j;^^±0.26 


BaBar [68] 


Semileptonic 


1.7±0.8±0.2 


BaBar 


Average [67] 


1.76 ±0.49 




Our average 


1.68±0.31 



Table 3: Experimental results for ^{B ^ t v) prior to this conference from Rosner and Stone [64]. 

was constructed assuming that the systematic uncertainties in the two measurements are uncorre- 
cted, while BaBar supplies their average. The measurements are all very consistent, perhaps too 
consistent, considering the quoted uncertainties. 

Prior to this conference one inconsistency stood out that could be the first sign of NR The CKM 
fitter [62] and UTfit [63] groups had explored the consistency of all processes, by relating the the 
CKM constants A, X, p and rj to individual measurements, and deriving values for these numbers. 
While the overall consistency was at the ~ 2a level, delving into the individual measurements 
showed a looming difference between the value of the CP violating asymmetry in — )• J /^rK^ 
decays, sin2j3, and the value of ^{B^ — T^v). The CKM fitter group showed this discrepancy 
nicely by predicting the values of sin2j3 and J^{B^ — )• T^v) without using either measurement 
and comparing the prediction with actual measurements; see Fig. 28. The experimental average for 
sin2j8 is consistent with the predicted value, while the pre-conference value of .'^{B^ — )• T^v) is 
higher by a bit more than 2a. This problem was taken quite seriously by some: one paper appeared 
entitled "Demise of CKM and its aftermath" [70]. As we shall see shortly, this conclusion may 
have been a bit premature. 

At this conference Belle presented a new measurement using hadronic tags [71]. In doing so, 
they made several significant advances. First of all they improved their low momentum tracking 
efficiency thus obtaining a factor of 2.2 improvement in the number of hadronic tags. Secondly, 
they increased the data sample by a factor of 1.7, they also changed their analysis technique ob- 
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Sin 2|3 



Figure 28: Allowed correlated values of 3§{B^ t^v) and sin 2/3, shown in the banana Uke shape, com- 
pared with the experimental averages, shown as the small black circle, prior to this conference. The pink x 
shows the newly measured value from Belle for hadronic tags only. In all cases the errors are a combination 
of statistical and systematic uncertainties. 




0.2 0.4 0.6 0.8 1 1.2 5 10 15 20 25 30 



Figure 29: Distribution of Eecl (left) combined for all the decays from Belle. The M^j^g, distribution 
(right) is shown for Eecl < 0.2 GeV. The solid circles with error bars are data, while the the solid histograms 
show the projections of the fits, with the dashed and dotted histograms indicating the signal and background 
components, respectively. 

taining an improved signal efficiency by a factor of 1.8, and using besides extra energy, Eecl> they 
include information from the missing mass squared, M^j^^,, calculated from the hadronic tag and 
the decay products. These distributions are shown in Fig. 29. 

The Eecl signal distribution peaks at zero but is widened by detector and initial state radiation 
effects. The M^j^^ signal distribution is easy to predict from well understood decays and is 
flatter than the background distribution adding to the discrimination power of the analysis. With 
these improvements the results is only 40% of its previous value: 

^{B- ^x-v) = (0.72;j^;2^±0.1l) x 10"^ , (5.2) 
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being only about 3a significant. This branching fraction is also shown in Fig. 28 and is completely 
consistent with the sin2j8 measurement. We await updates using semileptonic tags from Belle as 
well as BaBar updates. 

5.2 B^D(*'t-v 

The semileptonic decay B — >■ D(*)t^ v proceeds in the SM via the tree-level diagram shown in 
Fig. 5(a), where the c quark and light spectator quark form a D or D* meson. As in the previous 
case a new charged boson could also modify this diagram especially if its coupling to the were 
enhanced [72, 73]. 



MB^D** (r/r^p 
B Background 



"-2 2 4 6 8 

ML. (GeV^) 

Figure 30: BaBar data for M^^-^^ and \p*f, \ (see insets) for Mj^j^^ > 1 GeV^ for the D+i' and D*+£' channels. 
(Data exist but are not shown for the D^i^ and D*^i^ channels.) In the background component the region 
above the dashed line is from charge cross-feed, while the region below contains continuum and BB. 

BaBar has recently shown a new analysis of these decays [74]. This analysis also requires tag- 
ging, but here only fully hadronic decay tags are used. The is observed only in its fully leptonic 
decays to or e^, so three neutrinos are missing. The main variables used to discriminate signal 
from background are the M^j^j, calculated here using the tag plus the lepton, and the magnitude of 
the momentum of the lepton in the B rest frame, \p*^\. The data are shown in Fig. 30. The large 
shaded (blue) peaks correspond to either the B — )• D£ v or the B — D*l^v channel and thus in this 
case constitute a background. The solid black regions correspond to D** production and also are 
a background that has a similar shape to the signal. Only the non-shaded (red) region in the top 
plot and the non-shaded (green) region in the bottom plot constitute signal. While it is difficult to 
separate the D** from the signal components, BaBar has used events with an extra to estimate 
this contribution directly from the data. 

Results are reported in terms of ratios ^(D^*)) = 3S(B D(*)t"v) /,^(B D(*)^"v). These 
quantities are well predicted in the SM [73]. The comparison is shown in Table 5.2. Taken together. 
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SM theory 


BaBar value 


Difference 




0.297 ±0.017 


0.440 ±0.058 ±0.042 


+2.0a 




0.252 ±0.003 


0.332±0.024±0.018 


+2.7a 



Table 4: BaBar results for the ratios ^(Z)'*') compared to the SM prediction. 



SM theory Belle value Difference 
Table 5: Belle results for the ratios ^{D^*^) compared to the SM prediction. 



the results disagree with these expectations at the 3.4a level.^ 

Belle had previously measured these ratios [75] and found excesses beyond the SM predic- 
tions; their results are listed in Table 5.2. Note however that the Belle results are not published, even 
though they were disseminated in 2009, and have rather large uncertainties compared to BaBar. It 
would be really interesting to have an updated analysis using the new tagging techniques that Belle 
uses in their most recent B — )• z^v paper to see if NP in these channels can be estabhshed. 

6. Other searches: the dark sector and Majorana neutrinos 
6.1 The dark sector 

Could it be that there are 3 classes of matter? One being well known SM particles with 
charges [SU(3)xSU(2)xU(l)], the second being dark matter particles with "dark" charges, and the 
third being a form of matter that has both charges, called "mediators." 

In one such proposal [76] the mediator is a particle of the U(l) gauge group and is nicknamed 
a "dark photon," playing on the oxymoronic nature of the name. In one experimental study BaBar 
looked for such mediators coupling to b quarks [77]. The idea was to see if there were any events 
where the T(1S) decayed into a real photon plus no visible energy. No such events were found 
and limits on branching fraction between « 10^^ — 10^^ are set, depending on the kinematical 
distributions of the photon. 

The coupling between the dark sector and the quark sector is specified by e, and the dark 
photon mass is labeled as mA'. Several other searches have been made and reviewed by Echenard 
[78]. They are summarized in Fig. 31. There will be more experimental tests of these ideas, 
although these searches have already produced significant limits. 

^The two results are inconsistent the prediction of the Type II two Higgs doublet model, because they imply different 
values for the model parameters. 



0.297 ±0.017 0.701°: +2.0a 

0.297 ±0.017 0.48+°;22+0-°6 ^q^^ 

0.252 ±0.003 0.47+°;}^+°;°^ +1.8a 

0.252 ±0.003 0.48i[!}1+!!l!^ +1.8a 



24 



New physics from flavour 



Sheldon Stone 




10-^ 



-10 



10 



10-4 



10-2 



10-3 



1-5 



0.01 



0.1 

mA'(GeV) 



1 



10 



Figure 31: Constraints on the mixing parameters, £, as a function of the hidden photon mass derived from 
searches in Y(25,35) decays at BaBar (orange shading) and from other experiments [76, 79, 80] (gray 
shading). The red line shows the value of the coupling required to explain the discrepancy between the 
calculated and measured anomalous magnetic moment of the muon [81]. From Echenard [78]. 



6.2 Majorana neutrino production in B decays 

A crucial question in formulating theories beyond the SM is whether or not neutrinos are nor- 
mal Dirac spin 1/2 fermions, or if they are their own anti-particles as suggested by Majorana [82]. 
If neutrinos are indeed Majorana then they allow a process called neutiinoless double-j3 decay, 
shown in Fig. 32(a) [83]. This process can proceed for any value of the Majorana neutrino mass 
as the particle is virtual in this diagram. A similar process for B decays is shown in Fig. 32(b). 
Each jLt^ here can be replaced by an or a T^. Belle [84] and LHCb [88] have searched for these 



Figure 32: (a) Feynman diagram for nuclear double-)3 decay, (b) A related diagram for B 0^^*"^^ pL 
decays. 

decays without success. Upper limits are given in Table 6.2. The upper limit in the e^e^ mode 
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Mode Experiment Upper limit x 10 



B ^D+e e Belle < 2.6 

B-^D+e^lx- Belle < 1.8 

B-^D+n-jX- Belle < 1.0 

B ->D+ii LHCb <0.69 

B-^D*+^-^- LHCb <3.6 



Table 6: Upper limits at 90% CL on B decays to like-sign dilepton final states. 



is not competitive with nuclear double-j3 decay. The other modes, however, are unique since the 
measure coupling of the Majorana neutrino to muons. 

Majorana neutrinos with mass below that of the B can be searched for directly. Consider the 
annihilation diagram shown in Fig. 33. The process is similar to the one for B — V shown in 
Fig. 27, but here the neutrino being Majorana and on shell can transform into a /i^ and a virtual 
W+, which can materialize into hadrons. In particular the decay rate into 71+ and D+ has been 
predicted as a function of the coupling between the heavy Majorana neutrino and the lepton, |V^4| 
[86]. 5 




Figure 33: (Feynman diagram for nuclear double-/3 decay, (b) A related diagram for B /i 
decays. 

LHCb has searched for these decays. They did not find any signals but have set the best 
experimental limits on | V^4| as a function of the Majorana neutrino mass, for values just above the 
pion mass to just below the B mass, shown in Fig. 34 [88]. Other searches have been caiTied out 
using like-sign dileptons at higher masses [89]. 

7. Conclusions 

Although there is no compelling evidence yet for NP, Flavour Physics is very sensitive to po- 
tential effects at high mass scales. All NP theories must satisfy stringent experimental constraints. 
LHCb has been very effective at dispelling effects with marginal statistical significance, although a 

^Similar calculations have been performed for semi-leptonic decays [87], but these turn out to be less sensitive. 
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Majorana neutrino mass (MeV) 



Figure 34: Upper limits on |V'^4| at 95% CL as a function of the Majorana neutrino mass derived from the 
^ ji^ji^TZ^ channel. 

few remain. Will some be established when precision increases? Improving measurements in such 
decays B ^ K li^li^ , — > l^^l^^ , B — ^ T^v, B — > D^*it^v, CP violation in B^ decays may 
show NP effects with increasing precision, and need to be aggressively pursued. Flavour based ex- 
periments also are looking for evidence that neutrinos are Majorana, and direct links to dark matter. 
We are looking forward to defining the next theory beyond the SM either with our current and near 
term data or with new future facilities including Super-Belle [90], and the proposed LHCb upgrade 
[91]. 
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